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Abstract 

Tipping points associated with bifurcations (B-tipping) or induced by 
noise (N-tipping) are recognized mechanisms that may potentially lead 
to sudden climate change. We focus here a novel class of tipping points, 
where a sufficiently rapid change to an input or parameter of a system may 
cause the system to "tip" or move away from a branch of attractors. Such 
rate-dependent tipping, or R-tipping, need not be associated with either 
bifurcations or noise. We present an example of all three types of tipping 
in a simple global energy balance model of the climate system, illustrating 
the possibility of dangerous rates of change even in the absence of noise 
and of bifurcations in the underlying quasi-static system. 

Keywords: Rate-dependent tipping point, bifurcation, climate system 

1 Tipping points - not just bifurcations 

In the last few years, the idea of "tipping points" has caught the imagination in 
climate science with the possibility, also indicated by both palaeoclimatc data 
and global climate models, that the climate system may abruptly "tip" from 
one regime to another in a comparatively short time0 

This recent interest in tipping points is related to a long-standing question 
in climate science: to understand whether climate fluctuations and transitions 
between different "states" are due to external causes (such as variations in 
the insolation or orbital parameters of the Earth) or to internal mechanisms 
(such as oceanic and atmospheric feedbacks acting on different timescales). A 
famous example is Milankovich theory, according to which these transitions are 



^NB This version of the paper includes a correction to Section 2.1. 
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forced by an external cause, namely the periodic variations in the Earth's orbital 
parameters (see e.g. [H]). Remarkably, the evidence in favour of Milankovich 
theory still remains controversial, see e.g. f21| . 

Hasselmann [TO] was one of the first to tackle this question through simple 
climate models obtained as stochastically perturbed dynamical systems. He ar- 
gued that the climate system can be conceptually divided into a fast component 
(the "weather" , essentially corresponding to the evolution of the atmosphere) 
and a slow component (the "climate" , that is the ocean, cryosphere, land vege- 
tation, etc.). The "weather" would act as an essentially random forcing exciting 
the response of the slow "climate" . In this way, short-time scale phenomena, 
modelled as stochastic perturbations, can be thought of as driving long-term 
climate variations. This is what we refer to as noise-induced tipping. 

Sutera [23] studies noise-induced tipping in a simple global energy balance 
model previously derived by Fraedrich [S] . Sutera's results indicate a character- 
istic time of lO^yr for the the random transitions between the "warm" and the 
"cold" climate states, which matches well with the observed average value. One 
shortcoming is that this analysis leaves open the question as to the periodicity 
of such transitions indicated by the power spectral analysis [TU Fig. 7]. 

There is a considerable literature on noise-induced escape from attractors in 
stochastic models [3] ■ These have successfully been used for modelling changes in 
climate phenomena [22] , although authors do not always use the word "tipping" 
and other aspects have been examined. For example, Kondepudi et al jl3) 
consider the combined effect of noise and parameter changes on the related 
problem of "attractor selection" in a noisy system. 

More recently, bifurcation-driven tipping points or dynamic bifurcations [3] 
have been suggested as an important mechanism by which sudden changes in 
the behaviour of a system may come about. For example, Lenton et al [TOl [H] 
conceptualize this as an open system 

^=f{x,X{t)) (1) 

where X(t) is in general a time-varying input. In the case that A is constant, 
we refer to ([IJ as the parametrized system with parameter A, and to its stable 
solution as the quasi-static attractor. If X{t) passes through a bifurcation point of 
the parametrized system where a quasi-static attractor (such as an equilibrium 
point x{X)) loses stability, it is intuitively clear that a system may "tip" directly 
as a result of varying that parameter, though in certain circumstances the effect 
may be delayed because of well-documented slow passage through bifurcation 
effects [2]. Related to this, Dakos et al [4] have proposed that tipping points are 
recognizable and to some extent predictable. They propose a method to de-trend 
signals and then, examining the correlation of fluctuations in the detrended 
signal, they find a signature of bifurcation- induced tipping points. These papers 
have concentrated on systems where equilibrium solutions for the parametrised 
system lose stability, although recent work of Kuehn [15] considers tipping effects 
in general two timescale systems as occasions where there is a bifurcation of the 
fast dynamics. 
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The explanation of climate tipping as a phenomenon purely induced by bi- 
furcations has been called into question. For example, Ditlevsen and Johnsen [5] 
suggest that the predictive techniques to forecast a forthcoming tipping point [1] 
are not always reliable. Indeed, noise alone can drive a system to tipping with- 
out any bifurcation. Nonetheless, it seems that the ideas of bifurcation-induced 
tipping can give practically useful predictions, for example in detecting potential 
ecosystem population collapses [6] . 

In their review paper, Thompson and Sieber j26j discuss bifurcation- and 
noise-induced mechanisms for tipping. They examine stochastically forced sys- 
tems 

dx^ f{x,X{t))dt + g{x)dW (2) 

where W represents a Brownian motion. Using generic bifurcation theory, they 
distinguish between safe bifurcations (where an attracting state loses stability 
but is replaced by another "nearby" attractor), explosive bifurcations (where the 
attractor dynamics explores more of phase space but still returns to near the 
old attractor) and dangerous bifurcations (where the attractor dyamics after 
bifurcation are unrelated to what has gone before). In [25], Thompson and 
Sieber clarify that a timeseries analysis of bifurcation-induced tipping point 
near a quasi-static equilibrium (QSE) relies on a separation of timescales 

'tdrift ^ W^crit ^ '^stab (3) 

where Kdrift is the average drift rate of parameters, Kcrit is the decay rate for 
the slowest decaying mode of the QSE and Kstab are the remaining (faster de- 
caying) modes. However, it is not easy to define Kdrift in general (especially in 
a coordinate-independent manner) and there is no a priori reason for inequality 
(O to hold for a given system. 

Rate-dependent tipping has not previously been discussed in detail, but it 
has been identified in |28| as an important tipping mechanism that cannot be 
explained by previously proposed mechanisms {i.e. noise or bifurcations of a 
quasi-static attractor). This paper aims to better understand the phenomenon 
of rate-dependent tipping by introducing a linear model with a tipping radius 
and discussing three basic examples where this type of tipping appears. 

We suggest that tipping effects in open systems can be usefully split into 
three categories: 

• "B-tipping" where the output from an open system changes abruptly or 
qualitatively due to a bifurcation of a quasi-static attractor. 

• "N-tipping" where noisy fluctuations result in the system departing from 
a neighbourhood of a quasi-static attractor. 

• "R-tipping" where the system fails to track a continuously changing quasi- 
static attractor. 

We demonstrate that each mechanism on its own can produce a tipping re- 
sponse. Furthermore, any open system may exhibit tipping phenomena that 
result from a combination of several of the above. 
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Figure 1: The system {x,y) in (a) represents a (high dimensional) open system 
with inputs S,{t). We identify a low dimensional but nonlinear subsystem x 
forced by some "observables" X{y) from the high dimensional system. The be- 
haviour of X in (a) can be partially understood by examining the open subsystem 
(b) for a suitable class of temporal forcing A(i). 



The paper is organized as follows; in the remainder of this section we discuss 
a setting for open systems allowing discussion of the three types of tipping 
phenomena. In Section [5] we formulate a criterion for R-tipping. In Section [3] 
we discuss three illustrative low dimensional examples of R-tipping; two related 
to bifurcation normal forms and one for a slow-fast system. Section |4] gives an 
illustrative examples of all three types of tipping for an energy-balance model 
of the global climate for different parameter regimes. Section [5] concludes with 
a discussion and some open questions. 

1.1 Towards a general theory of tipping in open systems 

Dynamical systems theory has developed a wide-ranging corpus of results con- 
centrated on the behaviour of autonomous finite dimensional deterministic sys- 
tems - often called closed systems, because their future time trajectories depend 
only on the current state of the system. If the systems have inputs that can 
change the fate of system trajectories then we say the system is open. Real- 
world systems are never closed except to some degree of approximation, and a 
range of methods have been developed to cope with the fact that they are open: 
(a) One can view external perturbations as time variation of parameters that 
would be fixed for a closed system, (b) There are various theoretical approaches 
to stochasticity in systems, either intrinsic or external, (c) Control theory al- 
lows one to design inputs to control a system's outputs in a desired way, given 
(possibly imperfect) knowledge of the system. 

In Fig. [ija) we illustrate an arbitrary high dimensional system where we 
have identified a low dimensional subsystem that we wish to check for "tipping 
effects" . We do this by analysing the response of an open system ([T]) in Fig.[Hb) 
to possible time- varying inputs X{t). Fig. [5] shows some possible candidates for 



4 




Figure 2: Candidate inputs X{t) in the subsystem ([T]). These may include 
noise- and signal- hke components as in (a), or purely deterministic smooth in- 
puts/parameter variations as in (b,c). Tipping responses in the subsystem may 
occur in response to noise fluctuations (N-tipping), to passing through a bifur- 
cation point for the parametrized subsystem A = Ao (B-tipping) or as a result 
of too rapid variation (R-tipping) . 

the input A(t); we are interested in classifying those inputs that lead to a sudden 
change in x. This "tipping" may depend on details of the noise (N-tipping), may 
involve passing through a critical value of A(i) corresponding to a bifurcation 
(B-tipping) of the parametrized subsystem or may depend on the rate of change 
of X{t) along some path in parameter space (R-tipping). 

2 R-tipping: a linear model 

We use a simple model to explore R-tipping and to give sufficient conditions such 
that R-tipping does/docs not occur. Suppose that the system ([1]) for x S M" 
and parameter A has a quasi-static equilibrium (QSE) i(A) with a tipping radius 
R > 0. For some initial xq with |a;o — S:{X)\ < R we assume that the evolution 
of X with time is given by 

— = M{x-x{\)) for \x-x{\)\<R (4) 

where M is a fixed stable linear operator (i.e. |e^^*| as t — >■ oo). More 
generally, we consider a time varying parameter, A(t), that represents the input 
to the subsystem. If \x{t) — x{X{t))\ < R then we say that x{t) tracks (or 
adiabatically follows) the QSE x{X). If there is a to such that \x{to) — x{X{to))\ = 
R then we say the solution tips (adiabatic approximation fails) at to and regard 
the model as unphysical beyond this point in time. The tipping radius may be 
related to the basin of attraction boundary for the nonlinear problem ([T]), as is 
the case in Sec. 3(a,b), but it need not be, as is the case in Sec. 3(c) and [5S]. 
System (|4]) shows only R-tipping - because M is fixed there is no bifurcation 
in the system and no noise is present. Clearly the model can be generalised to 
include M and R that vary with X{t), and/or nonlinear terms. Equation (|4]) 
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can be solved with mitial condition a;(0) = xq to give 

x{t) = e*^*xo - /* e*^(*-^'Mi(A(s)) ds. 



If we assume that the solution is modelled by the linear system near the QSE 
for an arbitrarily long past and set u = t — s then the dependence on initial 
value decays to give 



x{t) ^ - / e'""M5(A(< - u)) du. (5) 

Ju=Q 

Assuming that M is invertiblc and exponentially stable (more precisely, we 
assume that |e^'^*A/^'''i;| as i — )■ for fc = 1, 2) and that the rate of motion 
of the QSE and parameter are bounded (more precisely, the derivatives d'x/dA' 
and d'A/di' for ^ = 1,2 are bounded in time) then ([5|) can be integrated by parts 
to give 



x{t) = - [e^^"i(A(i - u))] !° - e^"^(A(t - u)) du 

.In dt 



and so 

/■°° dr 

x{t) - i{\{t)) ^ - e^'^=-{X{t-u))du. (6) 
Jo 

Integrating again by parts gives 

x{t)-i{t) = M-^^{X{t ~ u)) ^ / e^''^Ar^^{X{t-u))du 

ClZ Jq UiL 

= L(t) + E{t) 
The linear instantaneous lag is 

L(i)=A/-iJ(A(0) (7) 

If we define the drift of the QSE to be the rate of change 

, ^ dx dx dX 

^^'^■-Tt=TxTt 

then the linear instantaneous lag is 

L{t) = M-^r{t). (9) 
The error to the linear instantaneous lag is 

1"°° d'^T 
E{t) = - / e^^"Af-i---4(A(i - u)) du. 
Jo dt'' 

which includes historical information. This can also be expressed as 

poo 

E{t) = e^^"Af-i [i"{X{t - u)){X'{t - u)f + £'{X{t - u))X"{t - u)] du. 

To summarize, the solution of (|4]) follows the QSE x{X{t)) with a linear instan- 
taneous lag term L(t) and a history-dependent term E{t). 
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2.1 A criterion for R-tipping with steady drift 

If ^ = r is constant in time then we say the system has steady drift and ^ 
simpUfies to x{t) — x{X{t)) = M~^r. In other words, one can verify that E{t) = 
and that 

\x{t)-i{X{t))\^\Ar'r\, (10) 

On writing the matrix norm ||-A/|| = sup^_^Q |Mw|/|u| we note that for any r ^ 
and invcrtible M we have 

\\M\\-^ ■ \r\ < \M-^r\ < \\M-^\\ ■ \r\. 

Wc can avoid R-tipping if \x{t) — x{X{t))\ = \M~^r\ < R and hence a sufficient 
condition on the rate of parameter variation to avoid R-tipping is that 

\\M-'\\-\r\<R (11) 

while a sufficient condition for R-tipping to occur in this model is that 

\\M\\-'-\r\>R. 

In the intermediate case, the path of parameter variation will determine whether 
or not there is any R-tipping. 



2.2 General criteria for R-tipping 

In the more general cas^ where r{t) varies with t we use ([5]) and the inequality 
|gA'/«^| ^ |jgMu|j|^,| gjyg |-]^g upper bound 

dx , 



\x{t)~i{t)\ < 



dt 



{X{t~u)) 



du. 



If we define 



,{t) = sup 

s<t 



sup 

s<t 



then this means that 

\x{t) - x{t)\ < rmax{t) 

Moreover, if M is stable then 

||e*^"|| <ce- 
for some real c, /3 > (see [12] ) , and so 



,Mm I 



du. 



\x{t) - X{t)\ < r,nax{t) -. 



(12) 



^Thc first part of this section is a revised and corrected version of the corresponding 
section that appeared in in Phil. Trans. R. Soc. A (2012) 370, 1166-1184, and is based on 
the pubhshcd correction prepared with the assistance of Clare Hobbs. 
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Hence, we can guarantee that @ avoids R-tipping by time t if 



c 



{t) < R. 



(13) 



max 



If M is scalar then we can choose c = 1, /3 = — and ([T5| reduces to 



On the other hand, if M is a matrix then we need a good choice of c and (3 in (|12p 
to make the tipping condition (jl3p optimal, but this depends on the matrix 
structure and not simply the norm; see for example the text by Hinrichsen and 
Pritchard [T^] and the elegant estimates of Godunov [371 Eq.(13)]. 

Note that the tipping radius approach does not seem to easily give any 
sufficient condition to undergo R-tipping in terms of rmaxif), comparable to 
sufficient condition for tipping with steady drift given the end of the previous 
section. 

One can define a natural timescale for the motion of the QSE as 



note that in general, combinations of dx/dt and dX/dt do not give timcscalcs in 
units s~^. For an R-tipping to occur, this natural timescale may be comparable 
to the slowest timescale (e.g. the reciprocal of the leading eigenvalue of M) of the 
parametrized system. The three examples in the next section have \dx/dX\ = 1 
and i? « 1 so wc expect R-tipping when \dX/dt\ w However, if \dx/dX\ « 

1/e, then clearly we can have R-tipping even when |dA/o?t| « e||A^ li- 
lt is possible to think of more general tipping problems by analogy with the 
"linear system and tipping radius" model discussed here. For example, for an 
open nonlinear system we consider an "effective tipping radius" that corresponds 
to how far the linearized system needs to be from a branch of QSE to ensure 
that the nonlinear system tips. There is however no exact analogy possible - the 
effective tipping radius may depend on the shape of the local basin of attraction, 
the nonlincaritics present and the exact path taken in parameter space. 

3 R-tipping: model examples 

We give three illustrative examples of R-tipping. The first two are based on 
normal forms for the two basic codimension one bifurcation that are generic for 
dissipative systems: the saddle-node and the Hopf bifurcation. The third is an 
example that uses a fast-slow system to show that R-tipping can occur even in 
cases where there is a single attractor that is globally asymptotically stable. 




(t) < R 



(14) 
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(a) < r < /i (b) r = /I (c) r > /( 



A 




Figure 3: Phase portraits of ([15]) for (a) < r < ^, (b) r ~ ^, and (c) r > /i, 
including the two quasi-static equihbria, Xa and is, and the two invariant hues, 
A and B. 



3.1 Saddle-node normal form with steady drift 

We consider the example system for x e M with parameter X{t) G M and drift r. 

(x + xf-^i, (15) 

(16) 

with fixed fi > 0. In the (x, A) phase plane of (fT5 |) -(fT6 l) . there are two dx/dt — 
isoclines given by Xa{ij) = {(x, A) G ffi^ : A = —^JTL — x} and Xs{y) = {(x, A) G 
: A = y/Jl— x}, that correspond to two QSE; a stable node and a saddle, 
respectively, for ([TC)) . Furthermore, if /i > r there are two invariant lines, one 
attracting 

A{^x, r) = {(x. A) G : A = -VJT^ - x}, 

and one repelling 

S(/i, r) = {(x. A) G : A = VJ^~T - x}, 

both with a constant slope dX/dx = —1 (Fig. [3]). The stability manifests itself as 
an exponential decay (growth) of a small perturbations about A{^, r) {B{fi, r)). 

If < r < /i then B{^,r) defines a tipping threshold. Initial conditions 
below B{iJ,,r) converge to A{fi,r) whereas initial conditions above B{ii,r) give 
rise to solutions x{t) — >• oo as < — >■ oo. If r = /i then A{^,r) and B{fj,,r) 
coalesce into a neutrally stable invariant line AB [Fig. I^Jb)] that disappears for 
r > fj. [Fig. [3jc)]. Hence, ioi r > fi there is no tipping threshold, meaning that 
trajectories for all initial conditions become unbounded as t — oo. 

Let us assume that the system is at (xo,Ao) at time t — 0. If the initial 
condition (xo, Ao) lies between A = — x and Xs(m) then the critical rate Tc is the 
value of r at which the r-dependent tipping threshold B{fi,r) crosses (xo,Ao). 



dx 

dX 

'dt 



9 



If the initial condition lies on or below the line A = —x then the critical rate rc 
is the value of r at which r) and A{^, r) meet and disappear. This gives a 
precise value for the critical rate as the following function of initial conditions: 

^ ^ f M - (Ao + xqY if -xq<Xq< -xq + y^, , 
\ H if Ao < -xo. 

We can approximate this result using the simple linear model (HI with the 
linearization at the QSE as M = -2^ so that \\M-'^\\ = \\M\\-'^ = l/(2y^). 
Clearly, the linear model with R = l^fji given by (basin boundary for Xa) 
overestimates because the linear attraction weakens on moving away from 
the stable QSE in the nonlinear problem. This can be overcome by choosing an 
effective tipping radius Rc- Comparing with (jlip . the system avoids tipping if 

|r| < 2^lRc 

which, for rc = fj,, suggests an effective tipping radius Rc = y^/^. Finally, 
owing to steady drift, this problem can be reduced to a saddle-node bifurcation 
at r = ^ in a co-moving coordinate system y = x + X. 



3.2 A subcritical Hopf normal form 

As a second example we consider 

f-n.-A), (18) 

where z = x + iy Cz C. For the subcritical Hopf normal form with frequency w 
we choose 

F{z) = {-l + iuj)z+\z\'^z. 

Note that the system ([18]) has only one QSE at z = A(t). Two cases of R-tipping 
that we consider are with steady drift ( these can be reduced to a bifurcation 
problem in another coordinate system) and with unsteady drift (where there is 
no straightforward simplification to a bifurcation problem). 



3.2.1 Hopf normal form with steady drift 

Consider (fTS]) with a uniform drift of the QSE along the real axis, at a rate r 
(which must be real): dX/dt = r. There is a critical rate rc at which the system 
moves away from the stable QSE. We can find this rc analytically by changing 
to the co-moving system for w ~ z — X, 

^=Fiw)-r, (19) 

where a stable equilibrium represents the ability to track the QSE in the original 
system. Setting w — |u)| e*^ and rewriting equation (|19p in terms d\w\/dt and 
dO/dt gives an equilibrium at (|we|,0e) that satisfies 

|u;e|^ - 2|u;e|^ + (w' + l)\we\^ - r^ = 0. (20) 
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In the (r, w) parameter plane, there is a saddle-node bifurcation curve (5* in 
Fig-SIa)) whose different branches are given by equation pO|) with 



.e|± = U^l±\/l-J(l+^2)j, (21) 

and join at cusp points at (r,a;) ~ (±(2/3)'^/^, ±(1/3)^/^) (not marked in 
Fig. |3^a)). Linearising about the stable equilibrium (|ii;e|,0e) reveals that the 
characteristic polynomial 

+ (2 - 4\w,f ) S + UJ^ + {\We\^ - iMWel^ - 1) = 0, 

has a pair of pure imaginary roots indicating a Hopf bifurcation when |i«eP = 
1/2 and > 1/4. In the (r, w) parameter plane, (disjoint) Hopf bifurcation 
curves (H) originate from Bogdanov-Takens bifurcation points {BT) at (r, w) = 
(±1/2, ±1/2), and are given by 



r 



2 



and LV^ > 1/4, (22) 



that follows from equation (|20| with jwep = 1/2. At BT, saddle-node bifurca- 
tion changes from super (solid) to subcritical (dashed). It turns out that the 
stable equilibrium for ([T^ . indicating the ability to track the QSE in the origi- 
nal system, disappears in a supercritical saddle-node bifurcation when < 1/4 
and becomes unstable in a subcritical Hopf bifurcation when > 1/4. Hence, 
for initial conditions within the basin boundary of this equilibrium, the critical 
rate is given by 



rc{io)^\ ±V KI'^-2h e|i + (c.^ + l)KP- if c.2<l/4, ^23) 

Again, we can approximate this result using the simple linear model Q with 
the linearization at the QSE as A/ = (-1, 5; -5, -1) so that p/^^H = = 
0.1961. Clearly, the linear model with a tipping radius i? = 1 given by the 
unstable periodic orbit (basin boundary for z) does not account for nonlinear 
attraction away from the QSE and for the spiraling shape of trajectories when 
7^ 0. Therefore, we choose an effective tipping radius Re- Comparing with 
(fTTj) . the system avoids tipping if 

|r| < 5.0990i?c 

which suggests an w-dependcnt effective tipping radius Rc{^j->) = rc(a;)/5.099. 

R-tipping that reduces to a bifurcation problem in a co-moving system should 
not be confused with B-tipping: observe that the bifurcation parameter r does 
not vary in time, and it is "the ability to track the QSE" , rather than the QSE 
itself, that bifurcates. 
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Figure 4: (a) Solid curves in a two-parameter tipping diagram for (|18p with 
steady drift indicate the critical rate rc{uj). The stable equilibrium for the co- 
moving system p^ . or the ability to track the QSE in the original system ((T5)) . 
(b) disappears in a saddle-node bifurcation or (c) destabilises in a subcritical 
Hopf bifurcation when r = rdi-u). 



3.2.2 Hopf normal form with unsteady drift 

We now consider ([T5| where we include a smooth shift of QSE between asymp- 
totically steady positions at z = to z = A, according to 



A), A(to) = A/2 



(24) 



where p > parametrizes the maximum rate of the shift, A > is the amplitude 
of the shift and <o is the time when the rate of change is largest. Integrating 
(HI) gives 

X(t) = A(tanh(A p{t - to) 12) + l)/2. (25) 



which implies the following parameter dependence on time: 
X{—t) — > 0, X{t) —J- A as t ^ oo and 



d\ dX A'' I 

'dt-Tt^^''~~ 



Near t = to this describes a smooth shift between the location of an asymptot- 
ically stable equilibrium from z = to z = A, and the maximum rate of the 
shift is A^yo/4 at i = to- Observe that there is no change in stability or basin 
size of the QSE as t changes. Fig. [5] shows typical trajectories starting at an 
arbitrary initial condition within the basin of attraction using fixed A and two 
values of p. As shown in the diagram, there is a critical value pc such that for 
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Figure 5: R-tipping for ((T^ ~ (P^ for A = 8 showing time evolution for (a) 
p — 4.76 and (b) p ~ 4.8 (recall that p scales the maximum rate of change) from 
an initial condition (x, y, A) = (0.4, 0.5, 0.0001). For p > p^ ^ 4.78 we find that 
system trajectories no longer follow the stable quasi-static equilibrium (shown 
by the dashed line) as they meet its basin boundary. 



p < Pc the system can track the QSE while for p > pc a tipping occurs near 

t = to. 

Jan Sieber [pers. comm.) has pointed out that this case may still be quantifi- 
able by numerical approximation of the pc that gives a heteroclinic connection 
from an (initial) saddle equilibrium at (z,A) = (0,0) to a saddle periodic orbit 
at {\z — A|, A) = (1, A) for the extended system (|18p and ([24|l . Such a connec- 
tion indicates pc for which the (initial) saddle equilibrium moves away from the 
basin boundary of the stable equilibrium at {z, A) = (A, A). 

3.3 A fast-slow system with i?-tipping 

A particularly interesting case of R-tipping can occur in slow-fast systems that 
have a (unique, globally stable) QSE near locally folded critical (slow) manifold, 
of which the recently studied compost-bomb instability is a representative [18[ 
I28j . Here, we consider a simple example 

a;2 + A-hxi(a;i - 1), (26) 

N 

-E^i' (27) 

n=l 

with odd A^, fast variable xi £ M, slow variable X2 € K, and small parameter 
< e ^ 1. A unique equilibrium for ([26)) - ((27)) . i(A) = (0, —A), is asymptotically 
stable for any fixed value of A, and globally asymptotically stable if A^ > 5. The 
slow dynamics is approximated by the one-dimensional critical (slow) manifold, 
S{X) = {{xi,X2) e : X2 = -A - xi{xi - 1)}, that has a fold, L{X) = 
(i,— A + j), tangent to the fast xi direction. If A^ > 5, the fold defines a 
tipping threshold that is not associated with any basin boundary. Here, S'(A) 
is partitioned into the attracting part, S'a(A) for xi < \, fold L{\) for xi = ^, 
and repelling part, Sr{\) for xi > i. 




dX2 



13 



3.3.1 The slow-fast system with steady drift 

Consider (PS )) " (P7|) with a uniform drift of the QSE, i(A(t)), in the negative X2 
direction at a constant rate 

-=.>0, (28) 

so that A becomes the second slow variable. There is a critical rate, Tc, at 
which ([26|) - ((28)) is destabilized, meaning that trajectories diverge away from 
the QSE for r > Tc. We can find this critical rate in the singular limit, e — 0, 
by setting e = in (|26p . differentiating the resulting algebraic equation with 
respect to i, and studying the projected reduced system [7]: 



dxi 
~dt 



N \ 

(2xi-l)-i, 

ri=l / 



(29) 



d\ , , 

Tt = ^' 

that aproximates the slow dynamics for p6)) - ((28l) on the two-dimensional critical 
manifold, S ~ {{x\^xi^X) ^ : X2 = —A — x\{x\ — 1)} (gray surface in 
Fig. [6]). Although ([29| -(|30 p is typically singular at the one-dimensional fold, 
L = {(xi,a;2,A) G M'^ : xi = i, Xi — —A + -j}, its phase portrait can be 
constructed by rescaling time 

dt r 

— = -(2x1-1) ^ t = - {2xi{s) -l)ds, 
dr Jo 

producing the phase portrait for the desingularised system |14) : 

N 



dxi 

dX 
d^ 



-r{2xi - 1), (32) 



and then reversing the direction of time on the repelling part of the critical 
manifold, Sr- In this way, we find that for < r < X]^=i ^ " trajectories for all 
initial conditions within Sa converge to a stable invariant line that is defined by 
a constant xi satisfying r = X^^^Li^ii raeaning that trajectories remain close 
to the QSE, X, for all time [r < Tc in Fig. [HJa) ]. However, for r > J2n=i ^ "i 
trajectories for all initial conditions within Sa reach the fold, L, where they 
"slip off" the critical manifold and diverge away from the QSE in the fast Xi 
direction [r > rc in Fig. [6l[a)]. Hence, system (p6)) -p8 | exhibits R-tipping and, 
for e ^ 0, the critical rate is 

N 



^2-". (33) 
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Xl Xl 

Figure 6: R-tipping in slow-fast systems with a unique quasi-stable equilibrium, 
X, and (gray surface) locally folded critical (slow) manifold, 5 = 5a U L U S'r , for 
(a) the steady drift problem ((26 |) ~ ([28| and (b) the unsteady drift problem ([26| - 
([27]) and where e = 0.01 and = 1. In (a), Eq. §^ gives = 5 and 

shown are trajectories for r = 0.4 < Tc and r = 0.6 > Tc. In (b), Eq. (|37p 
gives pc ~ 0.99 for the initial condition at the origin (black dot) and shown are 
trajectories for p = 0.7 < pc and p = 1 > pc- 



3.3.2 The slow-fast system with unsteady drift 

We now consider (P5)) -(P7 | with a nonuniform drift 

dX \ ,„,x 

that is a logarithmic growth, A(t) = In [p(t — to) + e'^^*"^] . where wc assume 
p > 0. Again, there is a critical rate, pc, at which system ([26 |) ~([27 | and (|34p is 
destabilized. The key difFcrcncc from the steady drift problem is that pc depends 
on the initial condition within Sa- This is because the desingularised system 



dxi 

dX 
dr 



N 

e-'p-J2x"i, (35) 

71=1 

-e-^p{2xi-l), (36) 



has a saddle equilibrium for all p > 0, corresponding to a folded-saddle singu- 
larity [SIM]: 



F ■ 



ixi,F,XFip)) ^ (^^,-\nf2^2-/p^ 



for the projected reduced system. One can use the theory developed in [28l 
Sec. 4] to approximate the critical value, pc- Given F, the eigenvector 



-q/p+^2+{q/p)^,iy 
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corresponding to the stable eigendirection of the saddle F for an 
mitial condition (cci^Oi^o) within Sa, and as far as e — > 0, the critical rate can 
be calculated using [551 Eq.(4.12)] to give 

Pc « P exp f Ao + --ll^^^l^^] , (37) 

\^ -q/p + + (q/p)^ J 

where p = q = ^n=i '^2~". Below the critical rate, the trajectory 

misses the fold, L, and approaches the QSE, x, as time tends to infinity [p < pc 
in Fig. m^b)]. Above the critical rate, the trajectory reaches L and diverges 
from the QSE in the fast xi direction [p > pc Fig. lEtb)]. Note that in this 
example, the critical rate of parameter variation is of the same order as the 
slow dynamics - only when the parameter variation is very slow with respect to 
the slow variable and there are three timescales is tracking guaranteed. In this 
sense, the rate dependent tipping occurs when the slow and very slow timescales 
are no longer separable. 



4 B-, N- and R-tipping examples in a simple 
climate model 

We present a simple climate model that independently show, under differing cir- 
cumstances, all three types of tipping. In its deterministic version, this is a "zero 
dimensional" global energy balance model originally introduced by Fraedrich [5] : 

c— = Ri-Rt. (38) 

The state variable T represents an average surface temperature of an ocean 
on a spherical planet subject to radiative heating. Eq. psp is a deterministic 
energy conservation law where the constant c represents the thermal capacity 
of a well-mixed ocean layer of depth 30m covering 70.8% of the Earth's surface. 
The incoming solar radiation R J, and outgoing radiation R 'f are modelled as 

Ri^^pIo{l-ap{T)), Rt=esA<jT\ 

Here Iq is the solar constant and the parameter p allows for variations in the 
planetary orbit, or in the solar constant. An ice-albedo feedback is introduced 
to link variations in temperature with changes of ice and thus of albedo ttpi 
Fraedrich [8] uses a quadratic relation 

ap{T) =02- b^T^, (39) 

where the parameters 02 > 1 and 62 control the albedo magnitude and slope of 
the albedo-temperature relation. The outgoing radiation term is obtained by 
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lo 


1366 Wm-^ 




1 


a 


5.6704 X 10-8 Wm-^K-'^ 


h2 


1.690 X 10-5 K-2 


c 


10^ kgKs--2 


0.1 


1.6927 


&SA 


0.62 







Table 1: Values of the constants and parameters for Eq. (|40p . 



the Stefan-Boltzmann law, where esA is the effective emissivity and a is the 
Stcfan-Boltzmann constant. With these choices (|55)) is written as |H1 Eq. 4.1]: 

^ ^ /(T) = c-ia(-T4 + &^t2 - rf^), (40) 

a = es^cr/c, 6^ = /^/o62/4es^(T, = -^/o(l - a2)/4esAa'. 

Table [T] shows the values of constants and parameters for the system at equi- 
librium. Sutera [23] reformulates Fraedrich's model to incorporate stochastic 
forcing: 

dT = f{T)dt + ^/l}dW, (41) 

with /(T) as in PO)) . where dW is a normalised Wiener (white noise) process 
such that {dW)'^ has dimension of time i, has dimension of 1/t and the variance 
of dW per unit time is v. 

For ^ larger than a critical value < /.Jc < 1, the deterministic system (|40p 
has two equilibria T+ (stable) and (unstable). A saddle- node bifurcation 
takes place at some fi — with < /ic < 1, where the two equilibria 
merge and disappear. Sutera [23] studies N-tipping in the stochastically forced 
system (j41[) for /x > /ic, as a function of the distance /i — /tc from the bifurcation 
value. Namely, they compute the exit time such that the process jumps over 
the "potential barrier" and falls irreversibly to 'ice-covered earth'. 

We illustrate in Fig. [7] three situations where the Sutera-Fraedrich model 
exhibits "pure" B-, N- and R-tipping independently; parameter values are de- 
tailed in Table [2] Panels (a-b) show an example of R-tipping, (c) of N-tipping 
and (d) of B-tipping. For case (a-b), we evolve the dimensionless parameter A 
according to the ODE dX/dt = pX{l - A) and set 62 = (1 - A)6^'"* + A6^"^' - 
for this figure we use initial values A = 10^ and T = 290K. The value of 02 is 
calculated to ensure that 6^ — 4c?p is held constant for the parameter groups 
defined in (|40)) . The constant p can be thought of as simply scaling the rate of 
passage from the initial to the final values given in the Table. 

5 Summary and conclusions 

It is of great practical importance to understand the theoretical mechanisms 
behind tipping phenomena in the climate system as well as other systems. We 
have proposed here that such mechanisms can be effectively divided into three 
distinct categories: bifurcation-induced, noise-induced and rate-dependent tip- 
ping, respectively denoted as B-, N- and R-tipping. In particular, we describe 
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Figure 7: Illustrations of trajectories for the Sutera-Fraedrich model (|4T|) show- 
ing the presence of all three tipping types for parameters in Table [2]- horizontal 
axis (years) vertical axis (Kelvin). The solid lines show system trajectories while 
the dash lines show the location of the QSE - the branch T+ is stable while 
the branch T_ is unstable in this model. Panels (a)-(b) show an R-tipping for 
a smooth change of parameters between two steady states. In (a) p = 0.18 the 
system returns to the QSE after a transient. In (b) p = 0.19 the system becomes 
unbounded indicating a critical value pc ~ 0.185 yr^^. (c) shows an example 
of N-tipping in the presence of noise of amplitude = 1.0 yr^^; (d) shows an 
example of B-tipping on decreasing /x uniformly from 1 at a constant rate. Note 
that in case (d) the two QSE coalesce at a saddle-node bifurcation. 
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Parameter (a-b) (c) (d) 

fi 1.0 1.0 decreases from 1.0 

at rate -0.0004 yr^^ 

62 (K-2) initial 1.690 x 10"^ 1.690 x 10"^ 1.04 x 10"^ 

final 1.8350 x 10"^ 

02 initial 1.6927 1.6927 1.2 
final 1.8168 

(yr-i) 1.0 



Table 2: Parameter values for simulations shown in Fig. [7l For case (a-b) 
we interpolate between the values given along the curve such that 6^ — Ad^ is 
constant at a rate proportional to p. In case (c) all parameters are fixed but 
noise is added, while in case (d) we impose a steady drift of the parameter n 
downwards. 



R-tipping, a mechanism that may be exhibited by (sub-systems of) the climate 
system independently of the presence or absence of the other types of tipping. 

In realistic models, tipping effects may be associated with a combination 
of the three mechanisms, and it will be a challenge to understand this more 
general case. For example B-tipping, usually associated with slow changes in 
a parameter, may turn into R-tipping upon increasing the rate of change for 
the parameters. However, as schematically illustrated in Fig. |8l completely new 
mechanisms may appear on increasing the rate, including the possibility that 
B-tipping may be suppressed for fast enough variation of parameters. Alterna- 
tively, the B-tipping may persist but an R-tipping mechanism may come into 
play before the B-tipping is reached. 

We emphasise that neither N-tipping and R-tipping require any change of 
stability. Hence there is no reason to assume that the techniques of [5] , based on 
a de-trended autoregressive model for B-tipping, should deliver useful predic- 
tions in such cases - as noted by [5], N-tipping is intrinsically unpredictable. We 
are investigating whether any novel predictive technique may be developed for 
R-tipping. Those cases of R-tipping that can be reduced to a local bifurcation 
in a co-moving system may be expected to be predictable using similar meth- 
ods; this includes the examples in Sec. 3(a), (b)(i) with steady drift. In more 
complex cases, pc may still be quantifiable by global heteroclinic bifurcations for 
an extended system, for example (fTSll and ([24|) or (fTB]) and ([24]) in Sec. 3(b)(ii). 

The classification proposed here may be applicable to a wide range of open 
systems under the influence of noise and/or parameter changes. Recent work 
of Nene and Zaikin [5D] suggests there may be interesting applications of rate- 
dependent bifurcation theory to determine cell fate. There are potentially many 
other application areas, from mechanics and ecology to economics and social 
sciences where tipping points are of interest. We suggest that this will be an 
area of significant mathematical development in the coming years. 
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Figure 8: Schematic illustration of different possible system behaviours on 
ramping the parameter A at differing rates (dashed arrows) through the re- 
gion A S [0, f]. (a) shows an example where there is a B-tipping for low rates 
(quasi-static) that disappears for high enough rates, (b) shows an example 
where there is no tipping for small enough rates but R-tipping for large enough 
rates, (c) shows both B- and R-tipping, but there is a range of rates where no 
instability appears. 
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